Random quantum magnets with long-range correlated disorder: 
Enhancement of critical and GrifRths-McCoy singularities 
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We study the effect of spatial correlations in the quenched disorder on random quantum magnets 
at and near a quantum critical point. In the random transverse field Ising systems disorder corre- 
lations that decay algebraically with an exponent p change the universality class of the transition 
for small enough p and the off-critical GrifBths-McCoy singularities are enhanced. We present exact 
results for Id utilizing a mapping to fractional Brownian motion and generalize the predictions for 
the critical exponents and the generalized dynamical exponent in the Griffiths phase to d > 2. 



The presence of quenched disorder is known to have 
a pronounced effect on quantum (or zero-temperature) 
phase transitions: new universality classes emerge, some 
of them with unconventional scaling properties, quantum 
activated dynamics and strong, so called Griffiths-McCoy 
(GM) singularities even away from the critical point. 
A well studied example is the random transverse Ising 
model, for which many results in one dimension exist 
[^-^. Many of the unusual properties in Id were re- 
cently shown to persist also in higher dimensions 
In this paper we investigate the the effect of long-range 
correlated disorder on the critical and off-critical singu- 
larities at quantum phase transitions. This issue has, to 
our knowledge, never been studied before, although it is 
known to have significant impact on various other phys- 
ical phenomena like classical (thermal) phase transitions 
[H , surface properties [p^ljlll ] , anomalous diffusion in dis- 
ordered environments |12(| and many other areas. 

In addition to the experimental relevance for the quan- 
tum Ising spin glasses an interesting relation of the 
random transverse Ising models to the non-Fermi liquid 
behavior of /-electron compounds has recently been re- 
vealed [Q. For this system the low temperature prop- 
erties of the interacting Kondo impurities have been 
mapped onto an effective quantum spin- 1/2 system with 
strong Ising anisotropy and random bond and transverse 
fields, which turn out to possess long range spatial corre- 
lations that decay like a power law with distance @,|l5| . 
This observation, in addition to the above mentioned gen- 
eral interest, motivates the study of the present paper, 
which investigates the effect of long range spatial corre- 
lations in the quenched disorder on the quantum phase 
transition as well as on the GM singularities. 

To be concrete we consider the Ising ferromagnet with 
transverse fields 



H 



(1) 



(> 0) and/or the transverse fields hi (> 0) are quenched 
random variables. The spins are located on the sites of 
a d-dimensional lattice and {ij) denote nearest neighbor 
pairs on this lattice (modeling short-range interactions). 

At zero temperature (T = 0) this model has, in any 
dimension d a quantum phase transition from a param- 
agnetic to a ferromagnetic phase at some critical value of 
the average ratio between bonds and fields [ln{hi/ Jij)]^,-^. 
Here and in the following [. . .]av denotes the disorder av- 
erage. The distance from this critical point is 5, such that 
for ^ > the system (|l]) is paramagnetic, for 5 < it is 
ferromagnetic. We can introduce local deviations from 
the critical point by introducing variables (5(r), which in 
Id are simply given by 5{r) = \n(hr/ Jr,r+i)- Thus one 
can discriminate regions in space that tend locally to be 
ferromagnetic (paramagnetic) even for 5 > Q {5 <Q). 

Here we wish to study systematically the effect of 
(isotropic) spatial correlations in the disorder that can 
be modeled with a disorder correlator G(r): 
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Uncorrelated disorder is described by choosing G(r) to be 
a delta-function. The Harris criterion for correlated dis- 
order |^,|l6| shows that any disorder correlator that falls 
off faster than r'^/'^ (i.e. G(r) ~ 0{r-P) with p > 2/^, 
where is the correlation length exponent for uncorre- 
lated disorder) does not change the universality class of 
the quantum critical point of model (|l]) with uncorrelated 
disorder. On the other hand for 



G(r) 



with 0<p<2/iy {< d) : 



(3) 



where ct^'^ are spin-1/2 operators and the interactions Ji^ 



where the last inequality holds generally for disordered 
system with uncorrelated disorder [|l^ (c.f. 10), the 
disorder correlations are relevant (and thus truly long- 
ranged) , the critical exponents become different from the 
uncorrelated case and the quantum critical point con- 
stitutes a new universality class, which we are going to 
explore in this paper. 
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First we consider the Id case for which we can derive 
most of our results in a rigorous way utihzing a mapping 
of the problem to fractional Brownian motion |l9LpO 
First we note that v — 2 for uncorrelated disorder^ 
which means that the long-range correlations are relevant 
for p < 1. The critical exponents and scaling relations 
can be determined by studying the finite size scaling be- 
havior of model (Q). As it is shown in ||] the gap /S.E 
(lowest excitation energy) of a chain of length L is given 

by Ai? « msrnsh]^Ylf=i f^i/ where and rris 
are the left and right surface magnetizations 
respectively. Thus, 



(4) 



Since S{i) are random variables with zero mean at crit- 
icality that are correlated according to (H) we con- 
clude that InAE scales like the transverse fluctuations 
of a correlated random walk. Since [{J2f=i H^)}^]av ~ 
L J-^ dr G{r) « L^^p we have for p < 1 



In AE oc L'^'-P^ with ^{p) = 1 - p/2 



Thus for long-range-correlated disorder the quantum ac- 
tivated dynamics scenario at the critical point is even 
enhanced and we get a new critical exponent ip. Note 
that for p > 1 one gets ip = 1/2, the result for uncor- 
related disorder |^,||. We checked this result by com- 
puting numerically the probability distribution Pl[AE) 
(see g for details), which we indeed found to scale like 
PL(lnA£;) - L-'^('''p(ln A£:/L'^('')) with ^(p) as in (|). 

The surface magnetization = {af) of a finite chain 
of length L (with the spin at site L being fixed) is given 

by m. = {l + Y.tlVu=l{h./J^,^+l?)r"^ §■ From 
this expression and (S) one sees that [In ms(i)]av ^ 
—XjPip)^ i.e. that the typical magnetization decays with 
a stretched exponential. Moreover, away from the criti- 
cal point {5 > 0) one has [lnTOs(L)]av ^ —L5 implying 
[ms{L,S)]typ oc exp(— L/^typ), where we defined typical 
correlation length that is seen to scale like ^typ ~ (5~'^'yp 
with ftyp — 1 independent of the correlation exponent p. 

On the other hand the average surface magnetization 
can be shown 21| to scale like the survival probabil- 
ity PsnrviL) of a random walk of L steps. This can 
be related to the first return time distribution Pf.r.t.(-^) 
of a fractional Brownian motion (with Hurst exponent 
H ^l-p/2e [1/2 Jl, c.f. (|,|)), which has been shown 
to scale like L"-^ [|2|. Since Psurv(i) = ■Pf.r.t.(i) 
L^~^ we get for p < 1 



,iL)U ^ L-^^'^P^ with Xsip)=p/2 
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For p> 1 one has the known result Xg = Ps/v = 1/2. 

From the analogy to fractional Brownian motion one 
can also derive the exponent v describing the divergence 
of the average correlation length when approaching the 
critical point, ^ ~ l*^! A non-vanishing distance from 
the critical point implies that the disorder configurations 



are such that they give rise to a non-vanishing average 
for the step width [5(i)]av — S, i.e. the fractional Brow- 
nian motion is biased. For d > (in the paramagnetic 
phase) the return time distribution has an exponential 
cut-off beyond a characteristic length scale f that scales 
like (51/^-1 H, which yields for p < 1 



U-S-''<^P^ with i^{p) = 2/p 



(7) 



From (||) and (0) one gets Ps {p) = 1 , independent of p. 
For p > 1 it is z/ = 2. The finite size scaling behavior 
of the average surface magnetization is then described 
by [mgjav = L~^^rh{L^/'^S). In Fig. 1 we show a corre- 
sponding scafing plot for p = 0.75. 
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FIG. 1. Finite size scaling plot of the surface magnetization 
according to the form [rrisjav = L~^^fh{L^^'^5) for p = 0.75 
using (^) and (j^. The data are averaged over 50000 sam- 
ples using a symmetric binary distribution for the couplings 
{Ji,i+i e {A, A"^}, here A = 5, hi = 1). 

Inset: The bulk magnetization exponent xt as a function 
of the value of disorder correlation exponent p estimated by 
evaluating numerically the average persistence exponent for 
the same type of correlated disorder in the Sinai-model (see 
text). The horizontal line is the value for uncorrelated disor- 
der Xb = {3- 75) /4 = 0.19098 . . ., the dashed line is p/2 and 
represents the asymptotic dependence of xi,{p) for p ^ 0. 



'-'^ such a chain is 



The bulk magnetization mt 
much harder to calculate, see [|5[|. The size dependence of 
the average bulk magnetization at criticality determines 
the last and remaining critical exponent Xb via 



[mb(L)]av ^ i-"^'") 



(8) 



In the case of uncorrelated disorder it was possible to pre- 
dict the exact bulk magnetization exponent Xh = (3i,/i> = 
(3 — \/5) / 4 (that will also hold for any p > 1) using a 
particular real space renormalization group [0, which, 
however, appears to be inappropriate for long-range cor- 
related disorder. In the limit p the difference be- 
tween bulk and surface magnetization will vanish due to 
the extreme correlation of the disorder. Hence we ex- 
pect Xfc(p) « Xs{p) = p/2 for p 0. To obtain the 
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full p-dependence we computed numerically the average 
magnetization for finite systems using the fermion repre- 
sentation 1^ and analyzed its finite size scaling behavior. 
We found only a slight decrease of Xb{p) with p, which 
presumably does not reflect the true asymptotic behav- 
ior that is much harder to reach numerically for corre- 
lated disorder than for uncorrelated disorder. Therefore 
we used the relation between the scaling properties of the 
bulk magnetization and the average persistence of a Sinai 
walker ||2l| , |2^ to compute the exponent Xb{p), as shown 
in the inset of Fig. 1. It confirms the asymptotic behavior 
Xb{p) ~ p/2 and the expected inequality Xb{p) < Xs{p) 
p3[ . In passing we note that the order parameter profiles 
m\r) = [{<Jr)lsN do not have the simple scaling properties 
reported in P4[ for the uncorrelated case. 

Away from the critical point the physical properties 
are controlled by strongly coupled clusters (i.e. segments 
that have locally a tendency to order ferromagnetically 
[^) giving rise to the so called GM singularities. In 
the paramagnetic phase {S > 0) the probability Pl{1) 
for a strongly coupled cluster of length I is propor- 
tional to Le~'/^ implying a typical length of such a 
cluster ^typ ~ ^InL, with ^ given by (0). With eq(||) 
for the gap we therefore get a typical energy scale of 
InAE oc <5/typ oc -S^~''''PhnL. Therefore we obtain for 
p<l and (5 < 1 

A£; - i"^' with z'{S,p) <x5^-'^/P (9) 

(and z'{p,5) — 25^^ for p > 1, (5 ^ 1). The generalized 
dynamical exponent z'{p,5) parameterizes all singulari- 
ties occurring in the GM phase P,p|,p5[ : e.g. the spin au- 
tocorrelation function at T = decays algebraically with 
z', Gioc(t) = [(af (T)cri(0))]av ^ T~i/^'; the local suscep- 
tibility diverges for T ^ when z' > 1, xioc ^ r^/^'"^; 
the specific heat has an algebraic singularity at T = 0, 
C ~ T^/^ • the magnetization in the presence of an 
external longitudinal field (in the x-direction) scales as 
M ^ H^/^ etc. We computed numerically ||^ the prob- 
ability distribution Pl{AE) which we confirmed to have 
a power law tail Pl{AE) oc AE-^/'^'+'^ with z' as in M). 

Concluding Id case we stress that (for fixed distance 6 
from the critical point) z' increases monotonically with 
decreasing S, i.e. stronger disorder correlations generate 
stronger GM singularities. Due to the nature of these 
singularities this tendency is a direct consequence of an 
increasing probability for large clusters for increasing dis- 
order correlations. 

Before we proceed to the higher dimensional case 
{d > 2) we describe briefly the infinite randomness dis- 
order fixed point (IRFP) scenario, originally developed 
for uncorrelated disorder [p7| , but one can generalize it 
for the present correlated case. This phenomenological 
theory involves three exponents: the lowest energy scale 
AE and the linear size L of a strongly coupled clus- 
ter are related via InAE' ~ L'^, its magnetic moment 
scales as ~ L'l''^ and its typical size when approaching 
the critical point, the correlation length, will diverge like 



^ ^ \5\ " . All bulk exponents can be expressed via ip, cj) 
and v, c.f. Xb — d — (f>ip, z/typ = 1^(1 — ip) and in the 
Griffiths phase z' cx 6'"'''. For the Id treated 
above, it is = 1/2, (f> = (V5 + l)/2 and v — 2 for 
uncorrelated disorder and p > 1, whereas for p < 1 we 
obtained ^ = 1 — p/2, = (1 — Xf,)/?/' (see Fig. 1) and 
u = 2/ p. The exponent relations are satisfied for the 
correlated and uncorrelated cases. 

In higher dimensions the IRFP scenario still holds 
p7|,pt, irrespective of the presence or absence of disor- 
der correlations. However, the exponents (^, (/>, v) will 
change. We can make precise statements on the change 
of these exponents for the case of random bond or site 
dilution, for which the quantum phase transition occurs 
at the percolation threshold p — pc (with p being the 
bond- or site concentration). In this case the physics 
is completely determined by the geometric properties of 
the percolating clusters Q, which means that {'ip,(f),i') 
can be expressed by the classical percolation exponents, 
which are the fractal dimension -Dperc of the percolat- 
ing cluster, the exponent /3perc determining the prob- 
ability for a site being in the percolating cluster, and 
the correlation length exponent t'perc, respectively. It is 
known |^,|l6| that the disorder correlations are relevant 
for p < 2/i^orc, in which case one gets z^(p) — Vpcrcip) — 
2/p, V'(p) — DpcTcip), which will increase with increas- 
ing correlations, since then clusters become more com- 
pact, and (l){p) = {d - Ppcrc{p)/vpcrcip))/Dpcrcip) whcrc 
/3pcic(p)/^'pcrc(p) is possibly independent of p In the 
GM phase this implies for the dynamical exponent 

^''-b-Pc(p)r'^--/^ (10) 

which increases with decreasing p, again confirming the 
general tendency that disorder correlations enhance the 
GM singularities. 

For the generic non-diluted case the exponents need 
not to be identical with the diluted case, however, one 
still has V — 21 p for p < 2/i/uncorr, according to a gen- 
eral argument given in Moreover, -0 increases with 
increasing disorder correlations, since its value is con- 
nected to the geometric compactness of strongly coupled 
clusters. Thus, the dynamical exponent z' ^ S~'^'^ grows, 
again enhancing the GM singularities. 

Regarding the recent experiments on /-electron sys- 
tems we would like to point out that there is evidence 
jisf that the spatial correlations in the metallic com- 
pound Ui_2:Thj:Pd3 decay like r~^. If we assume that 
for the transverse Ising model with uncorrelated disorder 
we have z/ — 2/d (as it is the case for d=l [^js) and d=2 
[Bl, and also for other random quantum critical points 
i 28| , p9| ) such a decay with p = 3 is the marginal case and 
instead of modifications of the above quoted exponents 
strong logarithmic corrections will appear. In a scaling 
theory for the marginal situation one has to replace L 
and ^ by LlnL and ^In^, respectively yielding for d—3 
and p=:3 for the gap, correlation length and dynamical 
exponent in the GM phase 
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(11) 



respectively, where "0 and are the critical exponents 
of the 3d system with uncorrelated disorder. Obviously 
these logarithmic corrections will make it very hard to 
extract the critical exponents i}} and v for instance from 
experimental data, and, furthermore, will apparently vary 
when approaching the critical point. 

Finally we would like to make a few remarks on related 
quantum magnetic systems: a) Quantum spin glasses ||2^ ] 
are expected to behave very similar as the random fer- 
romagnets with respect to the introduction of disorder 
correlations. The frustration caused by the random signs 
of the spin interactions is irrelevant at the IRFP |27j and 
therefore the universality class and the critical exponents 
are not changed, b) Random XY, e.g. in Id 



H 



'i+l 



(12) 



and XXZ or Heisenberg systems have different features 
in Id and in d> 2. In the latter higher dimensional case 
it seems that the quantum critical point is not an IRFP 
Ip?! , ^ , however, disorder correlations will certainly af- 
fect the critical properties. In Id we encounter the same 
scenario as for the transverse Ising systems, in particular 
for XY and XX chains |^ , since these are equivalent to 
two decoupled transverse Ising chains. Most remarkably 
the transverse and longitudinal correlations still decay 
with the same exponent (in contrast to the pure case), 
however, more slowly with correlated disorder. Moreover 
the term random singlet phase is now inapproriate when 
p < 1 , since then larger units than only pairs of spins will 
be strongly coupled. 

To summarize we have studied the effect of long-range 
correlations in the disorder on the quantum critical be- 
havior of random magnets. We have shown the relevant 
correlations generally enhance the critical and off-critical 
singularities, essentially because large strongly coupled 
clusters appear more frequently. For the random trans- 
verse Ising system in Id we reported exact values for the 
critical exponents for arbitrary disorder correlation expo- 
nent p, also for the diluted case in higher dimensions and 
argued how these results can be generalized to generic 
bond/field randomness in c? > 2. With respect to the re- 
cent experiments on f-electron systems we have pointed 
out the existence of strong logarithmic corrections that 
complicates the measurement of the critical and the dy- 
namical exponenent. Finally we generalized our results 
to XX and XY quantum spin systems. 
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